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Abstract. We calculate the fluxes of radio, hard X-rays and gamma ray emission 
from clusters of galaxies, in the context of a secondary electron model (SEM). In the 
SEM the radiating electrons are produced by the decay of charged pions in cosmic ray 
(CR) interactions with the intracluster medium, while gamma ray emission is mainly 
contributed by the decay of neutral pions. We specifically applied our calculations 
to the case of the Coma cluster, and found that the combined radio and hard X- 
ray fluxes can be explained in the SEM only if very small values of the intracluster 
magnetic field (B ~ O.lfiG) are assumed, which in turn imply a large energy density 
of the parent CRs. The consequent gamma-ray fluxes easily exceed the EGRET limit 
at 100 MeV. This conclusion can be avoided only if most of the hard X-ray emission 
from Coma is not produced by Inverse Compton Scattering (ICS). 



1. Introduction 

Clusters of galaxies have recently revealed themselves as sites of high energy processes, 
resulting in a multiwavclcngth emission which extends from the radio to the gamma 
rays and probably beyond. In this paper we refer to the Coma cluster, because of the 
wide evidences now accumulated for the presence of these non-thermal phenomena. 

The high energy processes which produce the observable radiations are due to 
the presence of a non thermal population of particles originating most likely from 
the cosmic ray sources in the cluster. The important role of cosmic ray electrons in 
Coma and in a few other clusters of galaxies is known since a long time because of 
the diffuse radio emission which extends over typical spatial scales of order ~ 1 Mpc. 
This radiation can be interpreted as synchrotron emission of relativistic electrons in 
the intracluster magnetic field. However, the combination of energy losses and diffusive 
propagation of these electrons makes their motion from the sources extremely difficult, 
so that it becomes a challenge to explain the spatial extent of the diffuse radio emission, 
unless electrons are continuously reaccelerated in the intracluster medium (ICM) (see 
|jj for a recent review). To solve this problem the SEM was first proposed in ^ ||: 
in this model CRs (protons) diffuse on large scales because their energy losses are 
negligible and can produce electrons in situ as secondary products of pp interactions 
with the production and decay of charged pions. 

The production of charged pions is always associated with the production of 
neutral pions which in turn result into gamma rays mainly with energies above ~ 100 
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MeV. The flux of gamma radiation and high energy neutrinos due to the cosmological 
distribution of clusters of galaxies was calculated in Q] and || , where the resulting 
diffuse neutrino background was also evaluated. Fluxes from single clusters were 
also compared with the upper limits on the gamma ray emission from the EGRET 
instrument onboard the CGRO satellite. 

More recently UV || and hard X-ray Q observations of the Coma cluster have led 
to the first detection of large fluxes at these wavelengths. Thier interpretation based 
on ICS of relativistic electrons off the photons of the microwave background radiation 
requires an intracluster magnetic field of B ~ 0.1/iG |?J. 

In this paper we calculated the multifrequency spectrum of the Coma cluster 
in the context of the SEM and compared our predictions with the results of recent 
observations. We find that also for secondary electrons, the radio and hard X-ray 
observations imply an intracluster magnetic field B ~ 0.1/xG and large energy densities 
in CRs. As a consequence, the flux of gamma rays above 100 MeV exceeds the EGRET 
upper limit [^J. We also discuss alternative scenarios which do not imply large CR 
energy requirements. 

The plan of the paper is as follows: in section 2 we describe the propagation of 
CRs in a cluster of galaxies; in section 3 we calculate the fluxes of secondary radio, X 
and gamma radiation and we present our conclusions with application to the case of 
the Coma cluster in section 4. 



2. Cosmic ray propagation in clusters of galaxies 

The propagation of CRs in clusters of galaxies was considered in previous works pi BL || 
where the effect of diffusive confinement of CRs was investigated. We summarize these 
results in the following. 

The CRs produced in a cluster of galaxies propagate diffusively in the intracluster 
magnetic field with a diffusion time, over a spatial scale R, that can be estimated as 
t ps R 2 /4D(E), where D(E) is the diffusion coefficient. Very little is known about 
the diffusion coefficient in clusters, but as argued in |Q for the bulk of CRs the 
diffusion time at large distances from the cluster center (R > R c , with R c ~ 1 Mpc 
the radius of the cluster) exceeds the age of the cluster for any reasonable choice of 
D(E). Assuming the following form of the diffusion coefficient, D(E) = DqE v . the 
maximum energy for which CRs are confined in the cluster is E c ps (R^/Doto) 1 ^ , 
where to is the age of the cluster, comparable with the age of universe. Despite of the 
large diffusion times in the cluster, CR protons do not suffer appreciable energy losses 
in the interesting energy range, so that the propagation can be simply described by 
a transport equation with no energy loss term (see, e.g., pof). Let us first consider 
the case of a single point source in the center of the cluster, injecting CRs with a rate 
Qp(E). As shown in |To[ the equilibrium number density of CRs with energy E which 
solves the transport equation can be written in the form 



D{E) 27H/ 2 r J r /r max (E) 



where r is the distance from the source and r max (E) = v/ 'AD(E)tQ. It is interesting 
to note that for r <C r max (E) the CR distribution tends to the well known time 
independent form: n p (E,r) = Q p (E) / (AnrD(E)) (see [|[ || for details). 

In the case of CRs injected homogeneously in the ICM, the equilibrium distribution 
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of CRs can be written as 

n p (E)=K^fp-y , (2) 

where V is the volume of the cluster, e tot is the total energy in CRs injected in 
the cluster and we assumed that the injection spectrum is a power law in the CR 
momentum p. The constant if is a normalization constant determined by energy 
conservation. Clearly this solution breaks down close to the boundary of the cluster 
and at high energies, where CRs are no longer confined in the cluster. 



3. Cosmic ray interactions and secondary electron emission 

The main interaction channel of CR protons in clusters of galaxies is represented by 
pp collisions with pion production. The decay of neutral pions produces gamma rays 
with energy above ~ 100 MeV, while the decay of charged pions results in electrons 
and neutrinos. The production of gamma rays and neutrinos from clusters of galaxies 
was recently investigated in j|, ||] . The secondary electrons produced by charged pions 
can play a fundamental role in the explanation of not thermal emission in clusters of 
galaxies, and here we describe this point in a greater detail. 

The 'primary electrons' models proposed as an explanation of radio halos in Coma- 
like clusters have serious problems due to the severe energy losses that make difficult 
the propagation of the relativistic electrons out to Mpc scales, where the diffuse 
radio halo emission is observed. This problem would be solved if electrons were 
produced or accelerated in situ. As initially proposed in [| this is the case for 
secondary electrons, generated in CR interactions with the thermal gas in the ICM. 
The same electrons would also produce X-rays by inverse compton scattering (ICS) 
off the photons of the microwave background. 

The calculation of the production spectrum of secondary electrons is explained 
in detail in p0| , and will be summarized here. For both radio and X-ray emission, 
the relevant electrons have energies above ~ 1 GeV, and it is important to have a 
good description of the pp interaction in a wide range of energies. For pp collisions 
at laboratory energy less than ~ 3 GeV, the pion production is well described by a 
isobar model, in which the pions are the result of the decay of the A resonance [l3| . 
For energies larger than ^7 — 10 GeV we use a scaling approach. In the latter the 
cross section for pp collisions depends only on the ratio of the pion energy E„ and the 
incident proton energy E p (x = E v /E p ) and can be written in the following form 

where fn(x) is the scaling function given in [|l(J for the case of charged and neutral 
pions. We refer to jl3| for a detailed expression of the cross section in the low energy 
case. 

The production electron spectrum at distance r from the cluster center, assumed 
to be spherically symmetric, can be easily calculated according to the expression 

( E e » r) = , * 2 n H (r)c • / dE, dE„ 



[ ' dE p F„ {E n ,E p )F e (E e , E„ , E n )n p (E p ,r) , (4) 

JE th (En) 
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where is the differential cross section for the production of a pion with energy E„ 
in a pp collision at energy E p (see for details), F e is the spectrum of electrons 
generated by the decay of a single muon with energy E p and n p is the spectrum of 
CRs. The function F e depends also on the pion energy because the muons produced 
in the pion decay are fully polarized, and this effect is taken into account here. The 
gas density at distance r is assumed to follow a King profile 



n H (r) = n 



-3/3/2 



(5) 



where, in the case of the Coma cluster we use uq = 3 x 10~ 3 cm -3 , ro = 400 kpc and 
(3 = 0.75 (we use here H = 60 km s _1 Mpc^ 1 ). 

The equilibrium electron distribution, n e (E e , r) , is achieved mainly due to energy 
losses, dominated at high energy by ICS and synchrotron emission and at low energy 
by Coulomb scattering. The effect of losses is to produce a steepening of the electron 
spectrum by one power in energy, at high energy, and a flattening by one power of 
energy, at low energy. In the next subsections we outline the calculations of non 
thermal radio, X and gamma ray emission from a cluster. 



3.1. The radio halo emission 

Electrons with energy E e in a magnetic field B radiate by synchrotron emission 
photons with typical frequency 

v = 3.7x 10 6 B fl E^Hz, (6) 

where is the value of the magnetic field in fiG. The emissivity at frequency v and 
at distance r from the cluster center can be easily estimated as 

M-Mft,r)(f) f, (7) 

\ / syn 

where (dE e /dt) syn denotes the rate of energy losses due to synchrotron emission and 
dE e /dv is obtained from eq. (||). The total fluence from the cluster is obtained by 
volume integration. 

Some simple comments can help understanding general features of the radio halo 
spectrum: for this purpose let us assume that CR are confined in the cluster and that 
the density of intracluster gas is spatially constant. If the injected CR spectrum is 
oc E~~< then the equilibrium CR spectrum is, within the distance r max {E) defined 

above, a power law oc E p ^ 7+r '- ) . Provided the electron energy is E e > 1 GeV, 
the production electron spectrum reproduces the parent CR spectrum, so that the 
equilibrium electron spectrum in the same energy region is oc E e ( 7+,,+1 ). From eq. (^) 
it is easy to show that r) oc z/~(t+''')/ 2 . The volume integration gives the observed 
spectrum: in the simple assumptions used here (complete confinement) the integration 
over the distance r at each frequency v must be limited by a maximum value 
Tmaxiy) oc so that the resulting radio halo spectrum is oc v^ 1 ! 2 , independent 

of the diffusion details. Repeating the same discussion, it is easy to show that the 
synchrotron radiation produced by CRs non confined in the cluster volume has a 
spectrum as steep as j/ _ (T+ r ')/ 2 . Therefore there is, in principle, a break frequency 
where the spectrum steepens from a power index 7/2 to a power index (7 + r/)/2. In 
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any realistic scenario, there is a smooth steepening of the spectrum. The presence of 
a density profile in the intracluster gas also contributes an additional small steepening 
of the spectrum. All these propagation effects are self-consistently considered in eq. 




3.2. Non thermal X-rays 

The peak energy where most of the photons are produced by ICS of electrons with 
energy E e is 

E x = 2.7E% (GeV)keV . (8) 

So, that the emissivity in the form of X-rays with energy Ex at distance r from the 
cluster center is: 

te^O-M*,,)^)^. (9 ) 

Here (dE e /dt)ics is the rate of energy losses due to ICS and dE e /dEx is calculated 
from eq. (|^). As for the radio halo emission, the observed fluence is determined by 
volume integration. 



3.3. The gamma ray emission 

As pointed out above, the production of neutral pions in pp collisions results in the 
generation of gamma rays with typical energy E 7 > 100 MeV. The emissivity of 
gamma rays with energy E 7 at distance r from the center of the cluster is given by 

A-E™ ax pE™ ax n (E r) 

j 7 (E 7 ,r) = 2n H (r)c / dE p F„ (E^E p ) pV *V 1/2 ,(10) 

where E™ m (E 7 ) — E 1 + iiv\ Q / '(AEj) . The function F n o is again calculated using the 
isobar model for proton energy less than 3 GeV and with the scaling model for energies 
larger than 7 GeV. The observed gamma ray flux is obtained by integration over the 
cluster volume. 

In (T^] the contribution of secondary electrons to the gamma ray flux through 
bremsstrahlung was also calculated. Since its contribution is small if compared with 
the contribution from pion decay, we neglect here the bremsstrahlung gamma ray 
emission. 



4. Application to the Coma cluster 

In this section we apply our predictions in the SEM framework to the case of the 
Coma cluster, for which multiwavclcngth observations are available. 

Two extreme scenarios of CR injection are considered here: a point-like CR source 
in the cluster center and a uniform CR injection distributed over the cluster volume. 
In both cases we assume that the CR injection spectrum is a power law in momentum 
with power index 2.1 < 7 < 2.4, covering the range of values expected for first order 
Fermi acceleration at shocks as well as for other CR acceleration mechanisms. 
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The diffusion of CRs in the cluster is described by a diffusion coefficient derived 
from a Kolmogorov spectrum of magnetic fluctuations in the cluster, according with 
the procedure described in || . In this case the diffusion coefficient can be written as 

2/3 

cm 2 /s (II) 

where l c is the size of the largest eddy in the Kolmogorov spectrum. 

The procedure used to evaluate the expected fluxes is the following: we first fit the 
spectrum of the radio halo as given in |l4| for 7 = 2.1 and 7 = 2.4. This allows to find 
the value of the absolute normalization of the injection CR spectrum in terms of an 
injection luminosity L p as a function of the average intracluster magnetic field B^. We 
carried out this calculation for B tl = 0.1,1 and 2. After calculating L p in this way, we 
determine the hard X-ray flux and the gamma ray flux above 100 MeV according to 
the expressions given in the previous section. Our results for the radio and hard X-ray 
emission for the case of a single source are shown in Figs. 1 and 2 respectively. The 
three panels refer to values = 0.1, 1, 2, as indicated. The normalization constants 
and the integral fluxes of gamma rays above 100 MeV in the same cases are reported 
in Table 1, where the gamma ray flux is compared with the upper limit obtained by 
the EGRET experiment || . The last column in Table 1 contains the ratio of the flux 
of gamma rays from electron bremsstrahlung compared with the flux of gamma rays 
from pion decay, as calculated in pfj| . 



D(E) = 2.3 x l0 29 E(GeV) 1/3 B~ 1/3 



20kpc 




Figure 1. Fluxes of radio radiation from the Coma cluster calculated in the SEM 
for 7 = 2.1 (solid lines) and 7 = 2.4 (dashed lines). The three panels are, from left 
to right, for B M = 0.1, 1, 2. 



Fig. 2 shows that a joint fit for the radio and hard X-ray fluxes is possible only 
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for low values of the magnetic field (B^ ~ O.f), which in turn imply large values of 
L p (see Table 1). As a consequence, the gamma ray fluxes easily exceed the EGRET 
upper limit from Coma. These results indicate that the SEM cannot explain the 
multiwavelength observations of the Coma cluster without violating the EGRET limit. 
This conclusion is not appreciably changed in the case of homogeneous injection of 
CRs. In this case, the value of 7 is fixed by radio observations and is 7 = 2.32. In 
order to fit the hard X-ray data at the same time, an ICM magnetic field ~ 0.1 
and a total energy in CRs ttot ~ 8 x 10 63 erg are required. This value, averaged over 
the age of the cluster, corresponds to a typical CR luminosity of L p sa 2 x 10 46 erg/s, 
which implies a gamma ray flux above 100 MeV in excess of the EGRET upper limit 
by a factor ~ 3. Therefore, also for a homogeneous CR injection, the SEM fails in 
explaining the radio and hard X-ray observations at the same time, without exceeding 
the EGRET limit. 

This conclusion has, however, further important consequences for other models 
too. In fact, we checked that the CR energy densities obtained in the present paper 
are comparable with those obtained assuming equipartition, as done in recent papers 
(e.g. pH ) to explain the radio, hard X-ray and UV observations. Therefore, the 
gamma ray limit applies to other models as well and forces the CR energy density 
in clusters to be some fraction of the equipartition value. Actually, as shown in ]15| , 
the present gamma ray observations put weak constraints on this fraction, but future 
gamma ray observations will definitely do better. We can envision at least two other 
arguments against equipartition of CRs and the ICM in clusters of galaxies: first of 
all, the most powerful CR sources typically present in clusters of galaxies (see |i| ||) 
allow to achieve a CR energy density equal to a small fraction (typically ~ 5%) of the 
equipartition value. Moreover, the magnetic field derived in Q] (B ~ 0.1/iG), which 
is the main reason to call for equipartition, is quite smaller than the equipartition 
magnetic field in the cluster, and it seems difficult to envision a scenario where CRs 
are in equipartition but not magnetic fields, in particular if the origin of CRs and 
magnetic fields in clusters are related each other. 



Table 1 . Summary of the fitting parameter values 
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So, at the present status of the debate, it is necessary to look for possible alternative 
explanations of the hard X-ray excess observed in Coma with the SAX satellite, since 
there is no stringent argument in favour of the ICS origin of such hard X-ray tail. A 
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Figure 2. Spectrum of the diffuse X-ray emission from Coma. The three panels 
refer to the same values of the IC magnetic field as in Fig. 1. The shaded area shows 
the best fit to the HEA01-A4 and GINGA thermal emission data (open triangles) 
at T = 8.21 ± 0.20 keV [18]. The OSSE upper limits [17] are indicated by the open 
circles. The SAX data [7] are indicated by filled squares. For this data set, errors are 
shown at 90% confidence level. Arrows and labels show, for each panel, the energy 
ranges in which the three different data sets are located. Predictions of the SEM for 
7 = 2.1 (solid lines) and 7 = 2.4 (dashed lines) are shown in each panel. 



possible alternative was proposed in |l6| where it was speculated that the presence 
of a 11011 thermal tail in the electron distribution could account for the X-ray excess 
through brcmsstrahlung emission. If this turns out to be the explanation of the hard 
X-ray excess, then the radio and hard X-ray fluxes become unrelated and magnetic 
fields of order > l^G, as the ones suggested by Faraday rotation measurements, would 
still be allowed. As a consequence, the CR energy density required to fit the radio and 
hard X-ray data would be of the same order of that predicted in || || . In this case, 
the SEM still remains a viable option for the origin of the Coma radio halo emission. 
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